High Field Specific Heat of 2D Quantum Spin System SrCu2(B03)2 
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We report measurements of tiie specific heat of the quantum spin liquid system SrCu2(B03)2 in continuous 
magnetic fields H of up to 33 T. The specific heat data, when combined with a finite temperature Lanczos diag- 
onalization of the Shastry-Sutherland Hamiltonian, indicates the presence of a nearest neighbor Dzyaloshinsky- 
Moriya (DM) interaction that violates the crystal symmetry for H — 0. Moreover, the same DM interaction is 
required to explain the observed ESR lines for H\\c. These results indicate that spin-lattice coupling needs to 
be included in any realistic description of this system. 

PACS numbers: 75.45.+j, 75.40.Cx, 05.30.Jp, 67.40.Db. 75.10.Jm 



SrCu2(B03)2 is a quasi-two dimensional spin system with a 
singlet dimer ground state yj. It is the only known realization 
of the Shastry-Sutherland model 101, and exhibits a sequence 
of magnetization plateaux at high magnetic fields f3,^. The 
unique behavior of this quantum spin liquid results from the 
interplay between two different fascinating aspects of strongly 
coiTelated spin systems: namely geometrical frustration and 
strong quantum fluctuations. The spin s — 1/2 Cu^+ ions 
that are responsible for the magnetism are grouped in dimers 
within planes of the tetragonal SrCu2(B03)2 unit cell, with re- 
spective intra-dimer or nearest neighbor (nn) and inter-dimer 
or next nearest neighbor (nnn) separations of 2.9 A and 5 . 1 A. 
The coupling constants are estimated to be J 80 K for nn 
and J' 50 K for nnn The geometrical frustration of 
the spin-lattice leads to very localized triplet excitations that 
have a tendency to crystallize at high magnetic fields. This oc- 
curs when the concentration of triplets reaches certain values 
that are commensurate with the underlying lattice, becoming 
incompressible upon formation of a gapped structure. The 
magnetization plateaux at Hpi = 27 T, Hp2 = 35 T and 
= 42 T are a direct consequence of spin superstructures 
forming at triplet concentrations 1/8, 1/4 and 1/3 respectively. 



Recent ESR experiments Oy, 121 1^ revealed the spin triplet 
excitation energy to decrease Unearly with increasing mag- 
netic field, extrapolating to zero at a field value (H — 22 T) 
very close to Hpi. However, on approaching Hpi, the experi- 
mental ESR data deviates from this linear extrapolation, indi- 
cating a level anti-crossing between the first triplet excitation 
and the ground state. The anti-crossing implies some mixing 
between two states with different magnetization along the 
tetragonal c-axis. This observation cannot be explained by the 
U{1) invariant models (which are symmetric under rotations 
around the c-axis) proposed in previous works, for which 
is a good quantum number. 

In this letter we argue that an intra-dimer Dzyaloshinsky- 
Moriya (DM) interaction, which violates the observed crystal 



symmetry at 100/^ < T < 395K f^, is required to explain 
the low temperature specific heat of SrCu2(B03)2 in magnetic 
fields _ff > 18 T. We also show that this interaction gives rise 
to the ESR transitions between the ground state and the single 
triplet excitations that are observed for H\\c [^J^^. These 
results suggest that there is a structural phase transition at low 
temepratures that lowers the crystal symmetry. Such a tran- 
sition could be driven by a strong spin-lattice interaction that 
would be a relevant ingredient to explain the magnetization 
plateaux which are observed in this system. 

The single crystal sample of SrCu2(B03)2 used in this 
study was grown by the floating zone technique. Stoichio- 
metric amounts of CuO, SrCOs and B2O3 were mixed, pre- 
annealed, and then annealed at 870 °C. Finally, the powder 
was regrinded, pelletized and annealed in O2 several times. 
Rods were formed by hydrostatic pressing and the growth was 
performed in a Crystal System Optical Furnace at a growth 
speed of 0.25mm/h in O2. No additional flux was applied 
lIOll. The measurements of the specific heat C{T,H) of 
SrCu2(B03)2 in continuous magnetic fields up to 33 T were 
performed on two oriented single crystal pieces of 12.34 mg 
and 13.92 mg. Both were measured with H applied along 
the tetragonal c-axis and within the ab planes. A calorimeter 
made of plastic materials and silicon was used, employing a 
thermal relaxation time technique optimized for rapid data ac- 
quisition iflUflill . The magnetization of a piece of sample 
of approximate dimensions 1.5 x 0.9 x 0.5 mm'^ was mea- 
sured as a function of field using a sample-extraction magne- 
tometer in a 400 ms, 45 T pulsed magnet provided by the Na- 
tional High Magnetic Field Laboratory at Los Alamos fl3ll . 
The small size of the sample, placed in good thermal con- 
tact with liquid '^He or ^He below T = 4 K, combined with 
the relatively slow field sweep of the magnet helped mini- 
mize magnetocaloric effects so as to achieve an isothermal 
experiment ll 1 4II1 . For characterization purposes, supplemen- 
tary Mz versus temperature T measurements and specific 
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FIG. 1: Magnetization vs field for SrCu2(B03)2 at different temper- 
atures between 0.6 K and 10 K, as indicated, reveal a gradual evo- 
lution of the magnetization plateaux. The square point on the 10 K 
curve was used to compare with SQUID magnetometer data in order 
to obtain the magnetization units. Inset A: Magnetic susceptibility 
measured at H = 4 T in a SQUID magnetometer (circles). The solid 
line is the susceptibility calculated with the FTL method. Inset B: 
Two copper dimers in the CuBOa plane where the coupling constants 
J (nn) and J' (nnn) are indicated. 



heat measurements were made at lower fields using a com- 
mercial Quantum Design™ MPMS (SQUID magnetometer). 
Meanwhile, numerical simulations of the Shastry-Sutherland 
model, with which the experimental data are compared, were 
performed on a 20-site square lattice using the finite tempera- 
ture Lanczos (FTL) method 

To describe the present system, we consider the following 
Heisenberg Hamiltonian on a Shastry-Sutherland lattice \2^: 

(iJ) (iJ)' 

+ 5] D.(SixSj)+ ^ D'.(SixSj). (1) 

Here, and (i, j)' indicate that i and j are nn and nnn 
respectively. The Hamiltonian includes nn (D) and nnn (D') 
DM interactions. The arrows indicates that the corresponding 
bonds have a particular orientation. The quantization axis z 
is parallel to the c-axis. The nnn DM interaction has already 
been considered in previous papers 1 18] to explain theposition 
of the single triplet excitations observed with ESR llala], far 
infra-red 11911 and inelastic neutron scattering measurements 
fl^. From the splitting between the two single-triplet exci- 
tations at q = 0, the D' = is estimated to be: D' =2.1 Kz. 
The orientation of the nnn bonds is given in Ref. d. Ac- 
cording to the crystal symmetry |9], the xy component of D is 
perpendicular to the corresponding dimer. However, as is ex- 
plained below, a non-zero z component of D, which is not al- 
lowed by the observed crystal symmetry at high temperatures 
(lOOi^T < T < 395i^), is required to explain the specific heat 
and the ESR data a a function of the applied field H. 

Figure n shows the magnetization Mz{H) measured as a 
function of magnetic field, in units of ji s/Cu determined upon 
cross-calibration with SQUID magnetometry data. In addi- 
tion to the plateaux already mentioned, there is a small ex- 
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FIG. 2: Measured specific heat divided by T (symbols) vs. T com- 
pared with the calculated one for the Hamiltonian of Eq. Q (solid 
lines) for the magnetic field along the c-axis: A) < < 18 T and 
B) 20 < il/ < 33 T. The parameters are the same ones used to com- 
pute the magnetic susceptibility. The dashed line is the calculated 
C/T for D = and = 22T. 



cess contribution to identified in our data over the en- 
tire field range. This additional source of magnetization has 
been observed before in SQUID magnetometry data |3], al- 
though for reasons which are unknown to us, it is absent in 
published pulsed magnetic field data. For our samples, both 
SQUID magnetometry data and low field pulsed magnetic 
field data evidence a finite excess susceptibility of approxi- 
mately 0.115 X 10^'' emu/mol Cu, probably due to crystalline 
defects. Good agreement with the expected magnetization 
values at the plateaux is obtained by subtracting this value. 
Better agreement is obtained by subtracting a scaled Brillouin 
function with an initial slope of0.14x 10^"^ emu/mol Cu and 
a characteristic temperature of 5 K (see Fig.Q. After either of 
these substractions, there remains a finite value of at very 
low temperatures and magnetic fields that increases linearly 
with H. 

In the inset of Fig. ^ we compare the measured magnetic 
susceptibility x{T) (after subtracting a small constant value of 
0.14 X 10^'^ emu/mol Cu) and the curve obtained with the FTL 
method that is described below. We get an excellent agree- 
ment for: J = 74K, J' = 0.62J,D=(2.2K,±2.2K,5.2K)(the 
sign is different for each dimer in the unit cell), D'=(0,0,2.2K). 
The values of the g-factors, g\\ = 2.15 and g±^ ~ 2.08, have 
been obtained from a comparison between our theoretical cal- 
culations |21J (see Fig. 13 and the ESR spectra |6, 8]. 

In Fig.|2j we show the specific heat divided by temperature, 
C(r, H)/T, for different values of the magnetic field applied 
along the c-axis. The primary feature in the low temperature 
specific heat is a broad anomaly centered at T = 8.5 K that 
is gradually depressed by increasing H. This anomaly has 
been attributed |22] to S'^ = dimer excitations. Here, how- 
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FIG. 3: Measured specific heat divided by T (symbols) vs. T 
compared with the calculated one for the Hamiltonian of Eq. Q 
(solid lines) for the magnetic field perpendicular to the c-axis:A) 
< H < 18 T and B) 22 < i^ < 33 T. The parameters are the 
same ones used to compute the magnetic susceptibility. The dashed 
line is the calculated C/T for D = and H = 25. 5T. 



ever, we observe a small shift in temperature as function of 
H, indicating the involvement of states with Sz ^ 0. For 
H > 12 T, a second anomaly develops at lower temperatures, 
which we attribute to single-triplet excitations situated 3 meV 
above the ground state in zero field. The Zeeman interaction 
causes this triplet states to move to lower energies with in- 
creasing H. Figure|3]shows similar results for H ± c. 

In Figs. 12] and |3j we also compare the experimental re- 
sults with the results of a numerical simulation of C(T, H) /T 
made using the FTL method 1 15, 16]. This method is based 
on the Lanczos procedure of exact diagonalization, and uses a 
random sampling over initial wave functions specially adapted 
for calculation of thermodynamic properties. All the results 
were computed on a tilted square lattice of = 20 sites. 
There are many advantages of this method over the conven- 
tional Quantum Monte Carlo (QMC) methods, which are as 
follows: first, the minus-sign problem that usually appears 
in QMC calculations of frustrated spin systems is absent; 
second, the method connects the high- and low-temperature 
regimes in a continuous fashion, enabling the entropy den- 
sity and specific heat (per unit cell) to be computed as ex- 
pectation values (i.e. s = kslnZ/N + {H)/NT, where 
Z is the statistical sum). The specific heat is then given by 
Cv = T{ds/dT) = kB{{H^) - {Hf )/NT^. The main lim- 
itation to the validity of the results originates from finite-size 
effects which occur when T < Tfs- The actual value of Tfs 
depends strongly on the particular physical properties of the 
system. For gapless systems, Tfs can be defined by way of 
the thermodynamic sum Z{T) = Tr exp(-(i? - Eo)/T), on 
condition that Z{Tfs) = Z* :$> 1 |[l6|]. In the present case, 
this condition can be relaxed (Z* > 1) owing to the existence 



of a gap in the excitation spectrum when J'/ J < 0.7 and 
to the almost localized nature of the lowest excited states - 
triplet excitations. By comparing results obtained on two dif- 
ferent systems with = 16 and = 20 sites, we estimate 
Tfs < 1 K. 

For H < 18 T, the agreement between theory and ex- 
periment is very good, regardless of the inclusion of the nn 
DM interaction. Finite size effects are also very small, for 
H < 18T, due to the localized nature of the single-triplet ex- 
citations |20|. When H approaches 22 T for H\\c (or 25 T for 
H _L c), however, the inclusion of this interaction is required 
to explain the measured C(T, H) /T at low temperatures. For 
H\\c, this is explained by the fact that D is the only interaction 
that violates the conservation of Af^, by mixing the = 
singlet ground state of _ff (D = 0) with the single-triplet ex- 
cited state with AI^ = ±1. This mixing becomes effective 
only when the energy difference between both levels is com- 
parable to |D|. For H ± c, the same type of mixing is pro- 
duced by the z-component of D. In other words, the level 
crossing that would occur in absence of the DM interactions 
is replaced by level anti-crossing. This can be seen in Figs.|2]3 
and|3j) where we also show the calculated C/T for D = 
and H = 22r (H = 25. ST) for H\\c (H ± c). In absence 
of the DM interactions, the level crossing generates a peak of 
C /T atT — which is not consistent with the experiment. 
In contrast, the level anti-crossing moves this peak to higher 
temperatures in agreement with the experimental data. The 
anti-crossing occurs for different values of H in the different 
field orientations due to the anisotropv of the y-factor and the 
fact that D is parallel to the c axis o^j. At high temperatures 
(T > 20 K), the specific heat data deviates from the theoreti- 
cal prediction owing to significant phonon contributions. 

For T - 10 K and fields H > 20T, there are small devi- 
ations between the experimental curves and the calculations, 
which can be attributed to the inter planar antiferromagnetic 
interaction J"/ J ^ 0.21 that becomes relevant when the con- 
centration of triplet excitations increases. 

The components of the nn DM interaction are constrained 
by the crystal symmetry at low temperatures |9]. According 
to this symmetry, the z-component of D must be zero. Al- 
though a non-zero DM vector of the form {d, ±d, 0) improves 
considerably the agreement between experiment and theory 
for H\\c, it does not reproduce our experimental data in the 
prox imity of Hpi= 27 T for _ff _L c. Ultrasonic experiments 
12411 indicate that in this region the lattice gets strongly dis- 
torted by the application of the magnetic field . We speculate 
that this lattice distortion increases the z-component of the 
DM vector to lower the magnetic energy (it increases the level 
anti-crossing). 

This speculation is further supported by the measured ESR 
spectrum for H\\c |^. We show the ESR spectrum as a 
function of H for = 5.2 (Fig.Hi) and = OK (Fig.gji) 
calculated with the Lanczos method 1*27]. More specifically, 
we are computing the dynamical susceptibility along the di- 
rection perpendicular to the applied field using the method 
introduced in Ref. |23]. As it is pointed out in Ref. 
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FIG. 4: Contour color plot for the ESR spectrum for H\\c calculated 
with the Lanczos method in a 20 sites cluster for a) Dz = 5.2 K and 
b) Dz = 0. The values of the other parameters are the same ones 
used to compute the magnetic susceptibility. The experimental data 
points in a) are from Cepas et al., [6]. 



constant |D| = 6.1 K, and a nnn Dzyaloshinsky-Moriya in- 
teraction constant |D'| = 2.2 K. A non-zero value of D^, that 
is not allowed by the observed crystal symmetry ?A H — Q, is 
required to explain both the specific heat data for iJ > 18 T 
and the observed ESR |6, 8] spectrum for H\\c. This sug- 
gests that a lattice distortion that lowers the crystal symmetry 
is induced at low temperatures. A more detailed comparison 
between the calculated ESR spectrum and the experiment will 
be presented elsewhere i2lll . 
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Magnetic Field Laboratory and Estonian Science Foundation 
grant No.4931. J.B. acknowledges the financial support of 
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netic Field Laboratory is supported by the National Science 
Foundation(DMR90- 16241), the Department of Energy and 
the State of Florida. 



the observed ESR trasitions between the ground state and the 
single-triplet excitations are not allowed by the observed crys- 
tal symmetry at iJ = @]. In Fig.|4t we show that these 
ESR transitions can be explained with a non-zero value of D^, 
while the corresponding ESR lines are not present if = 
(Fig.^J)). None of the other components of D or D' can re- 
produce these ESR lines. Based on these observations, we 
propose that the crystal symmetry is lowered at low temepra- 
tures due to a strong spin-lattice interaction. Since the lattice 
distortion depends on the applied field [E^], we expect to 
be an increasing function of H (although we used a constant 
value Dz — 5.2K for our calculation). 

With the exception of D, all other physical parameters 
used in the model were determined from previous experi- 
ments. The values of g\\ and g±^ are obtained from the ESR 
spectra f^, 'i']. By including D, however, we are able to ac- 
count simultaneously for the ESR spectra as a function of the 
applied field (see also Ref. |21]), the temperature dependence 
of the susceptibility, and the low temperature specific heat data 
for > 18 T 

In summary, we have measured the specific heat as a func- 
tion of temperature in continuous magnetic fields up to 33 T. 
An excellent fit to the C(T, H) /T data for both field orien- 
tations is obtained for a nn exchange constant J — 74 K, a 
ratio J' / J — 0.62, a nn Dzyaloshinsky-Moriya interaction 
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